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Abstract 

The quantization conditions which come from the requirement of 
the singlevaluedness of the odderon wave function are formulated and 
solved numerically in the 4 dimensional space of the odderon charge 
and the conformal weight h. It turns out that these conditions are 
fulfilled along one dimensional curves parametrized by a discrete set 
of values of Re /i in 3 dimensional subspace (Im /i, Im gs, Re (73). The 
odderon energy calculated along these curves corresponds in all cases 
to the intercept lower than 1. 

1 Introduction 



The leading contribution to the total elastic scattering amplitude of two 
hadrons (A, B) can be written in Quantum Chromodynamics (QCD) in a so 
called Regge limit 

s — > 00, t = const. (1.1) 

as a power series in a strong coupling constant of the partial amplitudes 
with a given number n of reggeized gluons (Reggeons) propagating in the t 
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channel: 



n=2 ^12) 

Here x„ denotes quantum numbers of the n Reggeon state and the residue 
functions Z?^'"^" and P'^'^ measure the overlap between the hadronic wave 
functions and the wave function of a compound state of n reggeized gluons. 
The n-Reggeon's partial amplitudes are proportional to s""(*) where a„(0) is 
called an intercept. 

It is of great importance to calculate the intercepts a„(0) in QCD, not 
only because they govern the high energy behavior of the forward elastic 
amplitudes but also because, e.g. for n = 2, they are responsible for the 
small Bjorken x behavior of the deep inelastic structure functions. The lowest 
non-trivial contribution for n = 2 was calculated in the leading logarithmic 
approximation by Balitsky, Fadin, Kureav and Lipatov 0121; who derived 
and solved equation for the Pomeron intercept. The equation for three and 
more Reggeons was formulated in Refs. [3 IH IH] • It took, however, almost 20 
years before the solution for n = 3 was obtained in Refs. [HI EI . 

The real progress started with an observation that the n-Reggeon ex- 
change is equivalent to an eigenvalue problem of a Schrodinger like equation 
with calculable interaction hamiltonian 7i„. Here the eigen-energy is related 
to Q;n(0) — 1. This Schrodinger problem is exactly solvable [11121110] which 
means that there exist n — 1 integrals of motion (g2, • • • , Qn) which commute 
with i-in and among themselves. The eigenvalue of q2 is equal —h{h — 1) 
where h is called a conformal weight. 

In the present work we shall concentrate the odderon exchange, i.e. on 
the case with n = 3. It is easier to conduct the calculations in the impact 
parameter space, that is in the transverse spatial coordinates of n Reggeons 
[xj, Uj). After introducing the complex coordinates {zj = xj+yj, z* = xj—yj) 
for j-th reggeized gluon, the odderon Hamiltonian becomes holomorphically 
separable 

n, = H, + % = ^Yl [^(^'^' ^^^+1) + ^('k, 4+i)] (1-3) 

k=l 

where is a number colors and zi = -Zn+i. The Hamiltonian 7^3 is con- 
formally invariant. Its eigenfunction is given as a bilinear form $ = ^ x ^ 
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where \I/ {^) is the solution of the Schrodinger equation in the holomorphic 
(antiholomorphic) sector. 

There have been many attempts either to directly find the values of £"3 
[TT]-[Tni or to find the spectrum of odderon charge gs [011121 • Finally, in Ref. [Zj 
the singlevaluedness conditions for the wave function $ were formulated and 
the spectrum of was found. This allowed to calculate the energy jB| (hence 
also the odderon intercept) for the conformal weight h = 1 — h = 1/2 which 
supposedly gives the largest contribution to the elastic amplitude. It is, how- 
ever, interesting to see explicitly whether h = 1 — h = 1/2 gives really the 
largest intercept and whether the there exist other solutions to the single- 
valuedness conditions than the ones found in [7j. The first attempt in this 
direction has been undertaken in Ref. J7j where the spectrum of ^3 for arbi- 
trary conformal weight but for non-physical quantization conditions h = h 
was found. In Ref. jlH] the spectrum of ^3 for the specific choice of the 
conformal weight h = 1 — h = 1/2 + iz/in the limit of small u has been 
studied. This result of Ref. |T2] was confirmed and extended for arbitrary u 
in Ref. (IT] . 

The values of the conformal the weight h and eigenvalues of ^3 form a 
four dimensional space. In the present work we construct an algorithm and 
numerical code which allows to find the points in the {h, q^) space which 
satisfy the physical quantization conditions h = 1 — h and singlevaluedness 
conditions of Refs.[71ll7j. It turns out that these points form one dimensional 
curves in the four dimensional space {h, q^). We have found families of curves 
which are numbered by discrete values of Re h = 1/2 + m/2 with m G 
3Z. Therefore for given Re h these curves are effectively embedded in 3 
dimensional subspace (Im h, Re ^3, Im q^). 

Applying the method from Ref. [H] which allows to calculate the odderon 
energy for arbitrary h and q^, we have calculated the odderon energy along 
the singlevaluedness lines qz{h). As expected, the odderon energy has a 
maximum for h = 1 — h = 1/2, and is always negative. Our numerical 
procedures are precise enough to find 17 values of ^3 for h = 1 — h = 1/2 
with 9 digits accuracy. 

The paper is organized as follows: in Section |21 following Ref. [T7j, we 
write the odderon equation in terms of variable ^ suggested in Ref.pH] and 
find its solutions around ^ = ±1 and 00. Recurrence relations for the for 
these solutions are collected in Appendix A. Next, in Sect. El we recapitulate 
the method of Ref.jH] and construct a singlevalued odderon wave function 
relegating the detailed form of the singlevaluedness constraints to Ap- 
pendix B. The resulting spectrum of q^ and h is calculated and discussed in 
Sect. 0] The numerical algorithm used in this Section is described in detail 
in Appendix C. Finally in Sect. El we calculate the odderon energy along 
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the singlevaluedness curves found in Sect. |3 Conclusions are presented in 
Sect. El 



2 Solution of the eigenequation for the odd- 
eron charge ^3 

2.1 Origin of the equation 

As already said in the Introduction it is possible to find a family of commuting 
operators qt which commute with the holomorphic n-Reggeon Hamiltonian 

Hn W 

[if„,gfc] = 0, A; = 2,...,n. (2.1) 

It follows that the Hamiltonian and the operators have the same set of 
eigenf unctions. In terms of the holomorphic coordinates qt have the following 
form: 

qk — ^ ^ ^ Zi\i2^i2i3 • • • ^ikh^h^h ■ ■ ■ ^ik^ l^-^) 

where k = 2, . . . ,n, Zjk = Zj — z^ and dj = d^y 

For the odderon case, n = 3, we have only 2 operators 



n>j>k>l 

q3= ^ -iZjkZkiZijdjdkdi. 

n>j>k>l>l 



(2.3) 



Following Ref.jH] we will use conformally covariant Ansatz for \P 

^{zuZ2,zs) = z''/^^{x), (2.4) 

where 

{Zi - Z2)iZi - Z3){Z2 - Zs) {Zi - Z2){zs - Zo) 

x = -, (2.5) 



{zi - Zoy{z2 - ^0)^(2:3 - 2;o)2' {zi - Zo){z3 - Z2)' 

h is a conformal weight and Zq represents an arbitrary reference point. A 
particular feature of this Ansatz ()2.4|1 is that q2 is automatically diagonal 



q2^{zi, Z2, Z3) = -h{h - l)^{zi, Z2, 2:3). (2.6) 
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In representation ()2.4|) the eigenvalue equation for ^3 takes the following form 



h 



[X 



2)(a; + l)(2x 



x{x — 1) 



-ip{x) 



+ 



h 



2x{x-l) - '-{h~l){x'^ 



X 



ip'{x) 



(2.7) 



+ 2x{x - l){2x - l)ip"{x) + x^{x - 1) = iq^^J^x). 



Equation ()2.7|) has been studied in Ref. [7j where the quantization condi- 
tions for g3 were found by introducing the singlevaluedness constraints on the 
whole wave function of odderon <P. The singlevaluedness conditions are much 
simpler when we rewrite equation ()2.7|1 in terms of a new variable suggested 
in Ref. I 



1 (x-2)(x+ l)(2x- 1) 



3^3 x{x 
Putting into we have |I7| 



Pkie - 1) 



{2.1 



viO = 0, 

(2.9) 



where 



{h + 2){h-3) 
6 ' 



ph 



h\h-3) 



27 



g3 

3^3' 



and qs is the eigenvalue of the operator ^3. 

As we shall shortly see the odderon equation ()2.9|) is less singular than 
Eq. ()2.7|) and the solutions of the indicial equation around ^ = ±1 do not 
depend on h. 



2.2 Solution of the odderon equation 

The odderon equation ()2.9p has three regular singular points at ^ = — 1, ^ = 1 
and ^ = 00. We shall solve this equation using the power series method. It 
is a third order ordinary differential equation therefore it has three linearly 
independent solutions. We can write them as a vector 



'«(^; qs) 



M^^qs) 
M^^qs) 
U3{^;q3) 



(2.10) 
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2.2.1 Solution of the equation around ^ = ±1 

Solutions of equation ()2.9|1 around ^ = ±1 have the following form [T7] : 



(±1) 



(2.11) 



n=0 



(2.12) 



where Sj are solutions of the indicial equation 

_ 2 _ 1 

Si ~ 2' "^^ ~ 3' '^3 — 

and do not depend on h. The coefficients u^"^^ are defined in Appendix A. 
2.2.2 Solutions around ^ = oo 

The solution of equation ()2.9p around ^ = oo has a more complicated form. 
In this case we perform a substitution ^ = 1/7] and then solve the problem 
around r] = 0: 



-^(1 - V0 + [2(1 - v')v - Ml - v'n 



+ 



du 
df] 



Ph- + q 
. V 



[2.13) 



u = 



The solutions of the indicial equation depend on the conformal weight h 



2h 



r2 



h 

1-3, rs 



h 
3' 



(2.14) 



and are identical as the solutions of the indicial equation in the case of 
equation ()2.7|) . Since r2 — rg is equal to an integer number, one of the 
solutions, u^^\^), contains a logarithm. The other two differences r2 — ri, 
ri — rs become integer as well if h is an integer itself. Therefore we have to 
distinguish several cases. 

For h ^ 7j and qs ^ the solution of ()2.9|1 in vicinity of ^ = 00 reads: 



u 



u 



S°°^(e;g3)=(i/0^'^E"ff(V0" 

n=0 

00 

^°°He;g3)=(i/o"E<^(ver 



(2.15) 



n=0 

00 



U 



(00) 



(e; gs) =(1/0^^ E 4^(1/0" + u^r\i; g3)Log(i/o 



n=0 
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where the logarithm Log{z) is defined as: 



Log{z) = In \z\ + zArg(z), |Arg(2;)| < n. 



(2.16) 



The coefficients entering ()2.15p are collected in Appendix A. 

The remaining cases, i.e. when qs = and/or G Z0, should be consid- 
ered separately. 



2.2.3 Solutions around ^ = oo for qs = 

It is easy to observe that in equation ()A.4j) for q^ = the term u 

tends to infinity. In this case the solution of equation ()2.9p should be 
constructed separately. With gs = and h ^ Z> the solution is given by 



(oo) 
3,0 



«^He;g3 = 0) = (i/e)^^E^5n^^^°^(V0'^ 



:2.i7) 



n=0 



The coefficients u. 



(oo;g3=0) 
i,2n 



are defined in Appendix A. 



2.3 Antiholomorphic sector 

For given h and q^ we can find the solutions of Eq. ()2.9p around all singular 
points by means of Eqs. ()2.11|) . ()2.15p and ()2.17|) . Analogously, we can con- 
struct the solutions in the antiholomorphic sector. Here instead of using the 
conformal weight h and charge q^ we use their antiholomorphic equivalents: 
h and Similarly to F^g. ^IJUif we write the three linearly independent 
solutions as a vector: 



■v2{C;q3) 



(2.18) 



3 Quantization conditions for the odderon charge 
3.1 Transition matrices 

Each of the solutions around ^ = ^1,2,3 where ^,2,3 = ±1, 00 has a convergence 
radius equal to the difference between two singular points: the point around 

^Solutions of the equation H2.13|l around ^ = cx) for /i e Z are not considered in this 
work. 

^ bar does not denote complex conjugation for which we us an asterisk. 
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which the solution is defined and the nearest of the remaining singular points. 
In order to define the global solution which is convergent in the entire complex 
plane we have to glue the solutions defined around different singular points. 
This can be done by expanding one set of solutions defined around in terms 
of the solutions defined around C,j for ^ belonging to the overlap region of 
the two solutions considered. Thus, in the overlap region we can define the 
transition matrices A, F, f2, where 



u 



■(oo) 



U 



U 



(-1) 
(+1) 



(3.1) 



{C;q3) = ^{q3)u^'^\C;q3)- 



Matrices Z\, F and f2 are constructed in terms of the ratios of certain 



determinants. 
Wrohskian 



For example to calculate the matrix F we construct the 



W 



-1) 
'(+1) 



u 



(e; qs) 
(e; qs) 



u 



u 



'(e;g3) 
(C; 93) 



(+1) 

3 

'(+1) 



U 



(e; 93) 
(e; qz) 
(e; 93) 



(3.2) 



Next, we construct determinants Wij which are obtained from W by replacing 
j-th column by the i-th solution around ^ = — 1, i.e. for i = 1 and j = 2 we 
have 



12 



ut'\i:q3) 



u 



u 



u 



'(+1) 



(e; 93) 
(e; qs) 



(-1) 

1 

/(-I) 



u 



//(-I) 



(C; qs) 
(e; q?,) 
(e; qs) 



u 



,(+1) 



(^; 93) 
(^; 93) 
(^; 93) 



(3.3) 



The matrix elements F^j are defined as 



(3.4) 



Matrix F does not depend on ^, but only on and /i. In a similar way 
we can get the matrices A and Q and their antiholomorphic equivalents: Z\, 
r, 72. 



3.2 Quantization conditions and singlevaluedness of 
the wave function 

The odderon charge ga is connected to its antiholomorphic equivalent by 

% = -qh (3-5) 
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where an asterisk over ^3 denotes complex conjugation. There exist two 
possible choices for q^: the one given by Eq. ()3.5p and a similar one with the 
plus sign. This follows from the fact that the eigenvalues of holomorphic and 
antiholomorphic Hamiltonian, £3 and £3, are symmetric functions of ^3 and 
^3 respectively |2Qj. Only one choice, namely leads to the nonvanishing 
solution of the quantization condition^ 

The odderon wave function can be written as |T7] 

^.7^,3^3 (^,^*) = ^'^H^*f^'v''iC-,q3)A{q,,q,m-,qs) (3.6) 

where ^ = C,{z). The wave function ^hEgsq^i^^ ^*) contains the solutions 
of equation ()2.9|) and its antiholomorphic counterpart, ^(C; ^3) and v{C,*;q^) 
respectively, and a 3 x 3 matrix ^(^3; ^3) "sewing" the solutions of the both 
sectors. 

The wave function ^ has to be singlevalued. This means that it should 
not depend on the choice of the Riemann sheet for the variables z and ^. In 
formula ()3.6|) the term uniquely defined only ii h/3 — h/3 E Z. 

This leads to the quantization condition for the conformal weight h 



h = -{n + m) + iu and h = -{fi — m) + iu, 



(3.7) 



where /i and v are real and m/3 G Z. The latter condition follows from 
the invariance under the Lorentz spin transformations. The normalization 
condition of the wave function requires that /i = 1 for the physical odderon 
solution. 

The fact that the wave function ^ should be singlevalued imposes certain 
conditions on the form of matrix A. It follows from Eq. ()3.6|) that for the 
solutions ()2.11|) around ^ = ±1 the matrix element A^i is multiplied by a 
factor (1 =p ^)**(1 =1= ^*)*^. This expression is singlevalued only if Si — sj G Z. 
For Si of Eq. ()2.12j) this is true only for i = j. Therefore the matrices A^^^'^ 
have a diagonal form 



a 
0/5 
7 



3. 13) 



a 





/3' 

i 







(3.8) 



For the solutions around = 00 (12.151 12.17j) . and for ^ Z, the matrix 
element Aji is multiplied by a factor {1/ S,Y^{1/ S,*Y^ ■ One should notice that 
solutions of the indicial equation r2, '"s and r3 differ by an intege: 



■^Note, that because of the factor i in the definition ^ (|2.8() . our ^3 has different sign 
than the one in Ref.j7|. 

and are solutions of the indicial equation in the antiholomorphic sector. 
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therefore terms which correspond to the elements A23 i do not vanish. 
Furthermore, terms with a logarithm appear in the solutions 1*3°°^ and ^3°°^ 
for gs ^ Z. One can see that when A23 = A32 then in the sum the ambiguous 
arguments of the logarithms cancel out. Moreover, for gs ^ Z the term 
which corresponds to the matrix element ^33 is not singlevalued. It contains 
a square of the logarithm which does not occur in any other terms. For this 
reason the element ^33(^3 7^ 0,^3 7^ 0) has to vanish. 

Thus the matrices A, defined around = 00, have the following form 



p 
Oar 
r 



A'^'^\q, = 0,qs = 0) 







a V 



(3.9) 



Substituting equation p.ip into the wave function ()3.6|) . one finds the 
following conditions for the matrices ^4(^3, ^3) 



A'iq,)A^°^\q,,qs)Aiq, 



r {li,)A^-'\q,,q,)r{q, 
lf[q,)A^+'\q,,q,)Q{q, 



: A(-)(g3,g3). 



(3.10) 
(3.11) 
(3.12) 



Each of the formulae 1)3.10113.1'^ is equivalent to a set of nine equations which 
can be conveniently written in terms of the following 4 vectors: 





a 




'a'' 




P 


a = 


/3 


, h = 




, c = 


a 
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i 




T 



d 



We can now rewrite equations ()3.1()f3.1ip3.12|l in the following form: 
• equation ()3.10|1 : 

— for g3 7^ as 



BupC= 0, Biou,c = 0, 



(3.13) 



for g3 = as 



B'upd 



Bdiagd — d- 



(3.14) 
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• equation ()3.11|) as 

— * 

CupS = 0, ClowO. = 0, CdiagO. = b, (3.15) 

• equation ()3.12p : 

— for 7^ as 

Dupb = 0, Dion,b = 0, Ddiagb = C, (3.16) 

— for 53 = as 

D'J=0, DiJ=0, D',,J=d. (3.17) 

Definitions of matrices B, C, D . . . are given in Appendix B. 

Equations ()3. 13113. l7|) have nonvanishing solutions only if the determi- 
nants of matrices with subscripts up and low are equal zero. Moreover there 
should exist the unique solutions of these equations: a, b, c, d which depend 
only on one free parameter which can be fixed by normalizing the wave func- 
tion As we shall show in the next Sections these two requirements fix 
uniquely the "boundary conditions" for the eigen equation of the operator ^3 
and allow to calculate its spectrum. 

4 Spectrum of the odderon charge ^3 
4.1 Eigenvalues of ^3 for h = 1/2 

Let us first discuss physical solutions found in Ref. jTj which correspond to the 
conformal weight h = 1/2 and Re ^3 = 0. In order to calculate the eigenvalues 
qs we have solved Eqs. (13.13113.171) requiring vanishing of the up and low matrix 
determinants. After that, we have also checked the uniqueness of obtained 
solutions. The results, also for unphysical values of ^3 with Im ^3 = are 
displayed in Table ^ Entries labelled from to 4, 12 and 13 agree with the 
ones of Refs. ^71 [7] , while the remaining ones are new. 

In fact the eigenvalues ^3 for h = 1/2 form a discrete set of points sym- 
metrically distributed on the real and complex axis in complex ^3 plane. 
Therefore in Table 14.11 only a half of the spectrum is displayed. It has been 
shown [HI [7j that only the imaginary values of ^3 are relevant for the odd- 
eron problem; real eigenvalues correspond to a wave function which is not 
totaly symmetric under the exchange of the neighboring Reggeons. There 
exists also one eigenvalue qs = which does not correspond to a normalizable 
solution (2^, see however |21j . 
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IN U. 


ya 


Nn 


HZ 


Nn 

IN U. 


ys 








6 


68.600522343i 


12 


1.475327424 


1 


0.205257506i 


7 


109.214406900i 


13 


12.947047037 


2 


2.343921063i 


8 


163.296192765i 


14 


44.413830163 


3 


8.326345902i 


9 


232.769867177i 


15 


105.872614615 


4 


20.080496894i 


10 


319.559416811Z 


16 


207.320706051 


5 


39.530550304i 


11 


425.588828106i 


17 


358.755426678 



Table 1: Eigenvalues of the odderon charge ga for h = 1/2 

4.2 Eigenvalues ^3 for the arbitrary conformal weights 

h 

As seen from Eq. fl3.7|) the conformal weight h depends on two parameters: 
m/?) & and z/ G M. In order to find the manifolds on which conditions 
()3.13ll3.l'7jl are satisfied we have extended the domain of h and allowed m to 
be any real number. This defines a four dimensional space of the conformal 
weight and the odderon charge [h, q^) G Mf^. In each group of equations 
()3. 13113. iTj) there are two equations which contain matrices with subscripts 
up and low. The determinants of these two matrices are complex, so we can 
define the following fourvalued functions: 

fs : (Re h, Im h, Re q^, Im gg) — > 

(Re(det Bup), Im(det Bup), Re(det Blow), Im(det Blow)) for qs 7^ 
(Re(det 5;,), Im(det K,), Re(det BiJ, Im(det BU) for q, = 0' 

(4.1) 

fc : (Re h, Im h, Re gs, Im q^) — ^ 

(Re(det C^p), Im(det Cup), Re(det Ciow), Im(det Ciow)), 

(4.2) 

fo : (Re h, Im h, Re gs, Im gg) — > 

(Re(det Dup), Im(det Dup), Re(det Diow), Im(det -D^ow,)) for gg 7^ 
(Re(det D;^), Im(det D'J, Re(det D;„^), Im(det D,',^)) for gg = ' 

(4.3) 
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Thus, in order to calculate the spectrum of the operator ^3, we should 
find common zeros of all functions /b, /c and fa: 



Furthermore one should verify the uniqueness of the solutions for ()3. 13113. l7j) . 



In Appendix C we have described the numerical algorithm constructed to 
find roots of Eqs. ()4.4|) . Our numerical findings can be summarized as follows: 

1. Although we have formally allowed m to be a continuous real parame- 
ter, the solutions of Eqs. ()4.4|) exist only for m/3 & Z, {e.g. Re h = 1/2 
or Re /i = 2). 

2. For the above discrete values of m, that is for fixed Re h, the solu- 
tions of Eqs. ()4.4j) form continuous curves in 3 dimensional subspace 
(Im h,Re g3,Im ^3). 

3. It turned out that each of 3 equations ()4.4j) yields the same set of curves, 
provided that the solutions of Eqs. ()3. 13113. l7jl are unique. 

In the following we discuss two sets of solutions to Eqs. ()4.4|l . namely for 
Reh = 1/2 and Re h = 2. 

4.2.1 Spectrum of for Re{h) = 1/2 

In Figure ^ we plot spectrum of the odderon charge ^3 as a function of h 
for Re h = 1/2. One can see two sets of curves: the ones in the plane of 
Re gs = and in the perpendicular plane of Im ^3 = and the line ^3 = 
which belongs to the both classes. For all these curves, except for gs = 0, the 
minimum of |g3| occurs for h = 1/2. Going away from this point the absolute 
value of qs increases monotonically. Minimal values of \qs\ correspond to the 
points listed in Table [T] We have not found any curve located outside of 
Re = or Im ^3 = planes. 

The curves in the plane of Re gs = have been earlier found in Ref . |22j 
and also in Refs. [d EH! • 

The spectrum of q^ has the following symmetries: 

1. Re gs -Re q^, 

2. Im g3 -Im g3, 

3. h^l-h 

which follow from the symmetries of the odderon equation ()2.9|) . Indeed, 
equation ()2.9|) is invariant under the transformation ^ — ^ and q^ —q^. 
The last symmetry follows from the exchange symmetry between the holo- 
morphic and antiholomorphic sectors. 



fB = 



fc = 



/d = 0. 



(4.4) 
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Figure 1: Odderon charge ga as a function of h for Re h = 1/2 Presented 
curves are described with point with the higher value of Igsj. 



4.2.2 Spectrum of for Re{h) = 2 

In Figure El we plot the odderon charge gs as a function of the conformal 
weight for Ke h = 2. Except of a line with gs = which is analogous to the 
case of Re = 1/2, the remaining curves have more complicated character. 
Still, they obey the following symmetries: 

2. Im h — Im h Re gs — > —Re gs, 

3. Im h — s> — Im h Im ga — — Im q^. 

These symmetries are similar to the case of /i = 1/2. 

The first one is connected to Bose symmetry — ^) and the others 

correspond to exchange of holomorphic and antyholomorphic sectors. For 
h = 2, that is, if the imaginary part of h vanishes, equation ()2.9j) has not been 
solved. In this case the eigenvalues gs have been obtained by an interpolation 
of the neighboring points which satisfied the quantization conditions. 
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Figure 2: Odderon charge gs as a function of h for He h = 2 



4.2.3 Case for = 

During the numerical computations we have noticed that similarly to the 
other matrices A, the matrix A^'^^qs = 0,q^ = 0) depends only on three 
parameters. The remaining ones vanish identically. In the case Re /i = 1/2, 
^' and v' vanish and for Ke h = 2, a' = t' = 0. 



5 Odderon energy 

The odderon energy is defined as 

Es = [ez{h, gs) + Qz)] = ^^^Re(e3(/i, gg)), (5.1) 

where £3 and £3 are the largest eigenvalues of the holomorphic and antiholo- 
morphic odderon Hamiltonian respectively. Applying the Bethe Ansatz we 
have for n = 3 
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for h=0.5 q3=0.205257i 
for h=0.5 q3=2.343921i 
for h=0.5 q3=1. 475327 
for h=0.5 q3=0.0 




Figure 3: Real part of holomorphic energy of odderon for Ke h = 1/2. Picture 
is plotted for curves from figure [T] 

where QsiX) satisfies the following Baxter equation 

(A + tfQ,{X + z) + (A - ifQsiX -i)- (2A3 + Xq, + q,)QsiX) = 0. (5.3) 
Equation ()5.3|1 was solved in Ref. [Hj by a substitution 



A*^Qfc(A) 



dz 
27x1 



K{z,X)P'Q{z)-zY, 



m=l " 



dz d 
2ni dz . 



z[z - l)L''-"'K{z, X)P"'-^Qiz) 
(5.4) 



where 

K{z,X)=z-'^-\z-iy^-\ 



L 



(5.5) 



Choosing the proper integration contour Cz and boundary conditions one 
arrives at a differential equation for Q{z): 

3 



.(.-1)^) -q.z\z-lf^^-^q,z{l-z) 



Qiz) = (5.6) 
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Similarly to the solutions of equation ()2.9p . the solutions of Eq. ()5.6p depend 
on the conformal weight h and the odderon charge ^3. Using the spectrum 
qzih) calculated in the last Section, we have calculated the energy of the 
odderon along the curves from Figures Q and El 

Analyzing the spectra of energy we can conclude that the odderon energy 
is always negative. This means that the intercept 

a{t = 0) = + 1 (5.7) 

is lower than one, so the odderon partial amplitude A^{s, t) ()1.2|1 is described 
by the convergent series in Regge limit (jl.lj) . 

5.1 Spectrum of the energy for Re h = 1/2 



53 


Re(£3) 




0.20526i 
2.34392i 
1.47533 


-0.73801 
-0.49434 
-5.16930 
-4.23462 



Table 2: Maximal values of Re(£:3) for Yie h = 1/2 

In Figure El we plot a real part of the holomorphic odderon energy Re(£:3) 
as a function of Im h. The picture is plotted for curves from Figure ^ that is 
for He h = 1/2. All these curves have a maximum in h = 1/2. The maximal 
values are displayed in Table El 

Going away from the maximum, the energy decreases monotonically. Our 
results agree with the values from Ref.j2|. The energy spectrum has the 
following symmetry 

h — >l-h. (5.8) 

5.2 Spectrum of the energy for Re h = 2 

In Figure^ similarly to Fig.jHl we plot Re(£:3) as a function of Im h for curves 
from Figure El that is for Re = 2. In this case the curves in {h, q^) space 
have much more complicated character and location of the energy maxima 
occurs not always for Im h = 0. 

For plots I, VI, Re(£:3) has a maximum ink = 2. In others cases the energy 
has a maximum in vicinity of h = 2, i.e. for the Il-nd curve the energy has a 
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-10,5- 

-1 1 ,0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-4-3-2-101234 

lm(h) 



Figure 4: Real part of the holomorphic odderon energy for He h = 2. The 
picture is plotted for the curves from figure |21 

maximum in {h = 2.0 + O.lOTi, = —3.508 + 2.050i). The maximal values 
of energy are given in Table El Similarly to the case of Re /i = 1/2 going 
away from the maximum, energy decreases monotonically. 

6 Summary and Conclusions 

The aim of the present paper was to look for the solutions of the odderon 
equation ()2.9p in the entire four dimensional space of the conformal weight 
h and odderon charge q^. So far these solutions have been found only for 
some specific values of /i [HI 123 ^| or for arbitrary h but unphysical 
quantization condition h = h ^Zj. To this end we have constructed and 
implemented the algorithm which is in detail described in Appendix C. This 
algorithm can be easily extended to more dimensional cases PE] . 

The calculations were performed in the holomorphic variable ^ which 
respects Bose symmetry and was proposed in Refs. |19| ITT] . The odderon 
equation ()2.9|1 is a third order ordinary differential equation with three regu- 
lar singular points at ^ = — 1, 1 and oo. Solutions around ^ = ±1 have been 
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curve symbols 


Re(£3) 


I 


-7.21 


II, III 


-6.27 


IV, V 


-8.51 


VI 


-4.00 



Table 3: Maximal value Re (£3) for Re /i = 2 
already found in Ref.^Tj. Here we have also calculated solutions around 

( = 00 dnsD- 

The singlevaluedness conditions imposed on the odderon wave function <P 
were found to be fulfilled along the discrete sets of continuous one dimensional 
curves. These sets are numbered by values ofRe/i = l/2 + m/2 (m/3GZ) 
and lie effectively in 3 dimensional subspace (Im h, Re q^, Im ^3). In this way 
we have obtained numerically the known spectrum of a Lorentz spin m. 

Although there are in principle 3 different singlevaluedness conditions 
obtained by gluing solutions around each of 3 singular points, which have 
to be fulfilled simultaneously, it turned out that it was enough to satisfy 
only one of them to get a complete set of solutions. It was therefore enough 
to consider two singular points, namely ±1 for which the solutions of the 
characteristic equation do not depend on h. 

Finally, we have calculated the odderon energy along the singlevaluedness 
curves qsih). For all cases the energy turned out to be negative which means 
that the odderon intercept is smaller than 1. The maximal value of the 
odderon energy corresponds, as earlier conjectured, to h = 1/2 and = 
0.205257506 x i which can be seen on Fig. 3. 

The authors thank J. Wosiek and G. Korchemsky for valuable comments 
and discussion. We are grateful to J. Wosiek and A. Rostworowski for making 
the program for calculating the odderon energy available to us. This work 
was partially supported by the Pohsh KBN Grant PB 2 P03B 019 17. 
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A Definition of series from solutions of the 
equation for ^3 

A.l Solutions around ^ = ±1 

The solutions of the equation ()2.9|1 around = ±1 have a form ()2.1H1 . where 
coefficients are defined by 



±1 1 ±1) ±1 / 

(A.l) 



,(±1) _i ,,(±1) _ . „,(±i) 
while 

aj,n-i = =F 4(ra + - 1) [{n + Sj - 2)(ra + Si) - /3J + 2ph ± 2g, 
hi,n-2 = -{n + Si-2) [{n + Si - 3)(n + s^) - 2(3h] + 2p/„ (A.2) 
mi,n =4(n + Si) [(n + Si){n + Si - I) + 2/9] . 

The upper sign corresponds to solutions around C = 1 ^"^^ the lower one to 
solutions around ^ = — 1. 

A.2 Solutions around ^ = oo 

The coefficients of solutions of equation ()2.15|) for gs 7^ and /i ^ Z we can 
write as 

• for i = 1,2: 



(00) (00) (00) / 



(00) M , / 



22 



u 



u 



u 



u 



u 



(oo 
3,0 

(oo 
3,2 

(oo 
3,3 

(oo 
3,4 

(oo 
3,n 



1 - h 



~ ' "3,1 



nlT^ = 0, 



2~q 

(«3,2?^3~^ + ^3,l4? + '^3,2^2^^ + f?,flU^2fi^ ^3,3, 

(«3,3M3~^ + h,2U^i^2 + C3,0M3^^ + d'i^^U^^t' + /3,lW2~^) /"^3, 

( (oo) I 7 (oo) I (oo) 

I '23,n-lW3 „_i + 03^n_2'U3^„_2 + C^,n-4:U^ n-4: 

+d^,n-lU'2^n-l + /3,n-3M£-3 + ^S.n-SM^-s) 



4, 



(A.4) 



where 



(A.5) 



0'i,n-i — 2g, 

&i,n-2 =2(n + Ti - 2) [(n + Ti - l)(n + r, - 2) - - 4/9] 

Ci,n-4 = - (n + - 4)(n + - 3)(n + - 2), 

c/i,„_i = - (n + r,) [3(n + r, - 2) + 4] + 2(1 + 

/,,„_3 =2 [(n + r, - 2)(3(n + r, - 4) + 8) - - 4/9] , 

5'j,n-5 = - 3(n + Ti - A){n + Ti - 2) - 2, 

+ Tifin + - 1) - 2(n + ri)(l + Ph) - 2ph. 



Similarly the coefficients for the solution ()2.17j) for qs = and /i ^ Z have a 
form 



(oo;<?3=0) 



(oo;g3=0) 
i,2n 



(oo;g3=0) _ , (oo) , 
1,2 — OiflUifl /mi,2, 

(oo;i33=0) 



l,oo;q3=u; . (oo;ij3=0) \ / 

'i,2(n-l)Mi_2(n-l) + Ci,2(n-2)^ii_2(n-2) ) /'^i,2n 



(A.6) 



B Definition of matrices connected to singl- 
evaluedness constraints on ^ 

In the section 1221 we defined quantization conditions for operator ^3. The 
matrices in formula ()3.13|) have a following form 



up 



Z\iiZ\i2 A21A22 + Z\21^32 

^11^13 ^21^23 ^31^23 + ^21^33 
^12^13 A22A23 A32A23 + Z\22^33 



(B.l) 



23 



low 



Bdi, 



diag 



Al^ll 
^13^13 



^21^21 
^22^22 
^23^23 



^31^21+421^31' 
^32 ^22 + ■^22 ^32 
^33^23 + ^23^33 



We can write matrices from ()3.14p as 



B' 



up 



B' 



low 





'AiiAi2 


A21A22 


^^31 ^22 


A21A32 


^31^32 




AuAi3 


A21A23 


^^31 ^23 


Z\21^33 


A31A33 


— 


^12^13 


A22A23 


A32A23 


zl22^33 


A32A33 






A22A21 


^Z2^21 


^22^31 


^32^31 






^23^21 


^^33^21 


^23^31 


^33^31. 




"AuAi2 


^21^22 


^31^22 


^21^32 


^31^32 




^11^13 


^21^23 


^31 ^23 


^21^33 


^31^33 




^13^12 


^23^22 


^33^22 


^23^32 


^33^32 




^12^11 


A22A21 


^32^21 


^22^31 


^32^31 




.^13^11 


^23^21 


^33^21 


^23^31 


^33^31. 




Al^ll 


^21^21 


^31^21 




^31^31" 




^12^12 


^22^22 


A32A22 


A22A32 


^32^32 




^^13^13 


^23^23 


^33^23 


^23^33 


^33^33 



B' 

diag 



The matrices occurring in ()3.15|) have a form 



up 



-^11 A2 

-^11^3 

BuBis 



-^21-^22 

-^21-^23 
-^22-^23 



-^31-^32 
-^31-^33 
-^32-^33 



Clow — Cup{r <— > r), 



c, 



diag 



-^12 A2 

-^13 A3 



-^21 Al 
-^22 A2 
-^23 A3 



-^31 Al 
-^32 A2 
-^33 A3 



Similarly we can write matrices from as 



D 



up 



'^^ll'i^l2 ''^21''^22 '^^31''^32 
i7lli7l3 i721^23 ^31^33 
^13^13 ^23^23 ^33^33 
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D 



low 



^12^11 



^22^21 ^32^31 

-^^23^^ -^^33^^ 

^23^22 — '^^22^23 ^33^32 — -^^32^33 



:B.ir 



diag 



^11^11 f22lf^2l '^^31^31 

^12^12 ^22^22 ^32^32 

(7213^212 + 7^12^?13)/2 (7223^222 + 7^22^?23)/2 (7233^232 + 7232^?33)/2 

(B.12) 

and the matrices from ()H.17|1 look like 



D 



up 



17111712 i^21^22 ^31^32 

17111713 i72ii723 ^31^33 
^12^11 ^22^21 ^32^31 



:b.i3) 



D' 



diag 



I22iJ721 ^31^31 

^12^12 ^22^22 ^32^32 

i7l3i7l2 i^23^22 ^33^32 

i7l2i7l3 i^22^23 ^32^33 

^13^13 ^23^23 ^33^33 



(B.14) 



:B.15) 



C Manifolds determined by set of M equa- 
tions in dimensional space 

C.l Method of finding zeros of the function F in N 
dimensions 



We shall be looking for solutions of a set equations 

Fj{xi,X2, . . . ,xn) = j = 1,2, . . . ,M. 



(C.l) 



Let X denote the vector of values Xj, and F the vector of functions Fj. Let 



us expand Fj{x) in a Taylor series 



N 



dF, 



Fj{xi + 8xi,xi + ... ,xn + Sxn) = Fj{xi, X2,. . . xn) + ^ 'g^^^'' ~^ ^('^^ 



i=l 



(C.2) 
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The matrix of partial derivatives appearing in equation ()C.2|) is the rectan- 
gular Jacobian matrix J. Thus, in matrix notation equation ()C.2|) reads: 

F{x + 6x) = F{x) + 36x + 0{6x'^). (C.3) 

We are interested in zeros of F{x), i.e. we are looking for such 6x that 
F{x + 6x) = 0. Neglecting terms of the order 0{6x we obtain a set of 
linear equations for the corrections 6x 

J6x = -F. (C.4) 

Equation ()('.4jl describes a set of M linear equations with n values of 
solution Sx. Each of these equations defines {N — 1) dimensional plane in 
dimensional space. In order to solve ()C.4|) we should find the intersection of 
these N — 1-planes. 

One should consider three cases: 

1. If the number of linearly independent equations is equal to the number 
of coordinates M = N then the set of equations has only one solution 

6x; 

2. If the number of linearly independent equations is lower than the num- 
ber of coordinates M < N and the equations are not contradictory, the 
set ()C.4|) has an infinite number of solutions which form a — M di- 
mensional plane. Then one selects the solution from this N — M plane 
which has the lowest norm. 

3. In the case when 6x doesn't exist, which means that the equations 
are contradictory, we adopt a procedure which tries to find some 6x' 
which decreases the test function (l(I3jl . If our set of equations is not 
contradictory then the algorithm reduces itself to two other cases. 

(a) Let m be a number of linearly independent equations which are 
not contradictory. Then, from the set of equations ()U.4|) we can 
construct k sets of m linearly independent equations. Each of 
these sets determines m entries of the dimensional vector Sx^^^ '. 
Here i (= 1, . . . ,k) corresponds to the i^^ set of equations. 

(b) The remaining (A^ — m) entries are set in a such way that Sx^^^' 
has the lowest norm. 

(c) Next we choose such j that sum of the angles between a vector 
6x^^^' and remaining vectors x^*-*' is minimal. 5x^^^' has the nearest 
direction to the average direction of other vectors x^*)'. 
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Instead of finding a zero of the M-dimensional function F we shall be 
looking for a global minimum of a function 

f = \F-F. (C.5) 

Of course there can be some local minima of Eq. ()C.5|) that are not solutions 
of Eq.jnil). 

Our step 5x is usually in the descent direction of / 

V/-(5f= (FJ) ■ (-J-^F) = -F-F< (C.6) 

which is true if only exists. 
It is convenient to define 

g{X)^f[XoM + \p) (C.7) 

where p = 6x. We always first try the full step i.e. X = 1. If the proposed 
step does not reduce / we backtrack along the same direction until we have 
an acceptable step 

Xnew = Xold + Ap, < A < 1, (C.8) 

i.e. we look for A which sufficiently reduces g{X). This is done by approxi- 
mating g{X) by polynomials in A. Initially we know g{0) and g'{0) and also 
g{l) which is known from the first trial (A = 1). Since we can easily calculate 
the derivative of g{X) 

g'{X) = V/ ■ p. (C.9) 
we can approximate g{X) by a quadratic polynomial in A: 

g{X) [g{l) - g{0) - g'iO)] X' + g'{0)X + g{0) (C.IO) 

and look for the minimum of ()(I10|1 . If this step also fails, we model ^'(A) as 
a cubic polynomial in A and so on until the satisfactory value of A is found. 
It is obvious that because of the linearity of the algorithm 

F{x + 6x) = 0{6x^). (C.ll) 

But, if f{x + 6x) < f{x), then our procedure leads towards the solution of 
()C.1|) . provided we are not in a vicinity of a false (i.e. local) minimum of /. 
In the latter case we have to change the initial conditions and start be whole 
procedure again. 
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C.2 Algorithm for finding the curves 



We shall describe a curve as a set of points placed along some path where the 
distance r between all neighboring points, should be constant. By definition 
all points should be zeros of function F. 

1. The input data of our algorithm are: 

- two points yi, y2, which should be placed in the vicinity of the 
sought curve, 

- distance r between these adjacent points. 

2. Making use of the algorithm from Section ITm we find a root Xi, situated 
in the vicinity of yi. 

3. We define the point xi as a center of a hyperspherical coordinate sys- 
tem. Next, we look on the sphere with radius r for the root X2 which 
has similar coordinates as point y2- 

4. We shift the center of the coordinate system to the point X2 and look for 
a zero of F on the hypersphere of radius r in vicinity of point yk{k = 3) 
extrapolated from previously found roots Xi with i = 1,2, . . . , k — 1. 

5. Iterating the above procedure we construct a curve of zeros of F. 

Using this algorithm we can find not only curves, but also k dimensional 
hypersurfaces. This can be done by fixing the values of — 1 coordinates. 
Then for the each choice depending on which coordinates were fixed, we 
can find curves from which we can in principle reconstruct the hypersurface. 
However, in our case it turned out that the zeros of F lie on one dimensional 
curves. 
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